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depends on the pulsar inclination which is clearly visible in the computed spectra of electrons 
escaping pulsar magnetosphere which become bimodal with the increasing inclination. This is an 
important result for modelling very high energy radiation from globular clusters where currently 
the standard power law or monoenergetic electron distributions are used. 
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1. Introduction 



Recent observations by the Fermi LAT have firmly established millisecond pulsars (MSPs) as a 
class of gamma-ray emitters. This long-awaited result is, however, accompanied by an astonishing 
finding: the radiation characteristics of all nine gamma-ray MSPs in terms of the phase-averaged 
spectra and the lightcurves are similar to the gamma-ray properties observed for classical pulsars. 
In particular, contrary to prior expectations, the results suggest that particle accelerators in MSPs 
are spatially extended in a way predicted by outer-gap and/or two-pole-caustic slot gap models. It 
is worth noting that currently all but one of the gamma-ray detected MSPs are pulsars from the 
field of the Galaxy. The exception is MSP B 1821-24, belonging to M28 globular cluster, which 
was detected by AGILE space observatory in 2009 ([|10|]). 

On the other hand, a multitude of MSPs reside in the dense cores of globular clusters (GCs). 
Due to a different evolutionary history in comparison with the field MSPs, one may expect that 
MSPs in GCs can have different structure of the accelerating gap than the field ones. For this reason 
we investigate, following [^] and [12], a millisecond pulsar radiation model where accelerating 
electric field extends from the polar cap surface up to very high altitudes in the volume filled with 
open magnetic field lines. The large spatial extent of the electric field is justified by a low number 
of created pairs, which is insufficient to effectively screen the electric potential. Such scenario is 
called a pair starved polar cap model (PSPC). 

On the basis of recent Fermi results ([JlJ], [j|]), it has been suggested that cumulative high energy 
emission of MSPs found in globular clusters contributes to their overall gamma-ray radiation. In 
addition to photons, MSPs inject relativistic electrons into GCs environment. These electrons 
propagate in the cluster and interact with its magnetic field producing X-ray radiation. Moreover, 
they are expected to up-scatter ambient photon fields (cosmic microwave background, infrared 
background and starlight) to high (HE) and very high energies (VHE). Recent theoretical studies 
show that such ICS photons may substantially contribute to HE and VHE radiation from globular 
clusters ([§, [JTTJ] , [§). 

In our study we investigate both the characteristics of HE and VHE radiation from MSPs and 
the distribution of relativistic electrons that escape pulsar magnetosphere. The results are based on 
a small sample of modelled MSPs, nonetheless, they show interesting characteristics of the ejected 
electron spectra which should be taken into account in modelling of the ICS component in GCs. 



2. Numerical simulations 



Calculations were performed using 3D numerical model of pulsar magnetosphere described in 
the framework of space charge limited flow. The magnetic field is in the form of a retarded vacuum 
dipole ([Q]) for which the curvature radii of magnetic field lines are determined in the inertial frame 
of reference. Moreover, following [Q] the special relativity effects like aberration and time-of-flight 
delays are treated accordingly throughout the calculations. 

An important assumption is that acceleration of particles takes place in the whole volume 
determined by the open magnetic field lines (PSPC model, e.g. [^]). The electrons are injected at 
the polar cap surface at the Goldreich- Julian rate with low initial energy. The accelerating electric 
field is in the form proposed by ^ which takes into account the effect of dragging of inertial 
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frames. In order to describe the electric field that extends up to the light cylinder three different 
formulae determining E\ | at different heights above the pulsar surface are combined. At moderate 
heights and close to the light cylinder the electric field is described by 



£ far ~E 2 exp 



+ E 3 (2.1) 



where E 2 is given by equation (14) and Ej by equation (35) of [^], T] = r//? NS , r is a radial distance 
from the neutron star and /?ns is a neutron star radius. In this approach T] c is a radial parameter 
which is determined for each magnetic field line via matching procedure. Its value depends on 
pulsar parameters like a magnetic field strength at the polar cap B pc and an inclination angle a (i.e. 
the angle between the magnetic axis and the rotation axis) and also on a position of the magnetic 
field line's foot-point at the polar cap. In the neutron star vicinity E neaT , given by equation (12) of 
[g], describes the electric field. We require the electric field Eu to be negative for all r up to the 
light cylinder. If this condition is not fulfilled for a certain magnetic field line, the line is left out 
from the calculations so no acceleration of particles takes place along this line. 

Our treatment of the PS PC model is different from that of [|l2|] in the way the near electric 



field ^near is matched with the far electric field Ef W . [ |12| ] first transit from E nsm to E% at a specific r\ 
equal to 175 « 1 -|__p-°- 333 . Then E 2 is matched with £3 so for distances close to the light cylinder 
the electric field converges to £3. To determine the electric field across the pulsar magnetoshpere 
we do not use the formula for T]b- In our calculations the matching of E new with Ef w takes place 
where condition £near — £far is fulfilled. This results in a different electric field structure across the 



pulsar magnetosphere with respect to [12]. 



3. Results 



The values of the key parameter (T7 C ) determining the shape of the accelerating electric field 
along open magnetic field lines are presented in the left panel of Figure [j] as a function of the line's 
foot-point at the polar cap. It is clearly visible that with the increasing pulsar inclination angle the 
structure of the electric field changes drastically. For small inclination angles most of the volume 
determined by the open magnetic field lines is filled with the electric field solutions described by 
high values of r\ c . As the pulsar inclination becomes higher the electric field solutions determined 
using low values of rj c start to appear on the part of the polar cap that is further away from the 
pulsar rotation axis {left panel of Fig. [I]). For large a (e.g. a = 70°) there is a part of the open 
magnetosphere where electrons are not accelerated because no T] c is found to match the electric 
field formulae. 

Such behaviour of the electric field structure with inclination angle has an influence on the 
spectrum of electrons escaping the millisecond pulsar magnetosphere. In the right panel of Figure |] 
different electron spectra are presented. For MSP with low inclination the electron spectrum is 
single peaked and the energies of electrons contributing to the peak are in the range from 10 6 to 
10 6 ' 7 MeV. The electron spectra of highly inclined MSPs are bimodal. As can be seen from Figure |T] 
electrons that are accelerated along lines with high value of the parameter T] c contribute to a high 
energy component (the same as for the low inclination MSPs) while those whre T] c is low build up 
a low energy component that streaches down to energies of 10 5 MeV. 
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The electric field structure affects also the photon characteristics of MSPs. Figure || shows 
photon maps and the light curves for selected observers in the energy range between lOOMeV and 
30GeV. Within the framework of PSPC model one obtains single pulse lightcurves. In the photon 
map for a pulsar with the inclination of 50° a wing-like feature is present. This corresponds to 
emission from the magnetic field lines anchored within the polar cap notch. 

4. Conclusions 

Our studies give a new insight into possible spectra of relativistic electrons that are ejected 
from the millisecond pulsar magnetospheres into globular cluster. The spectra are bimodal provided 
a > 20°, which is a new result with respect to the standard power law or monoenergetic electron 
distributions that are currently used in modelling of the VHE emission of globular clusters. 
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Figure 1: Presented plots are for MSP with P = 5 ms, B pc = 3.5 x 10 G and the inclination angle a, going 
from top to bottom: 20°, 50° and 70°. Left panel: schematic view of electron distribution across pulsar polar 
cap. Colour coded are different values of a parameter rj c used for combining the near electric field solution 
with the far solution (green = 5.6, red = 1.1, blue = all values between 1.1 and 5.6). Small and big solid 
circle mark the distance across the polar cap equal to 0.5r pc and 1.0r pc , respectively. Here r pc is the polar 
cap radius. Horizontal red dashed line shows <J> pc = 0° while vertical one points to <t> pc = 90° (Ope is the 
magnetic azimuth angle). Right panel: spectra of electrons escaping MSP magnetosphere. Colour coding is 
the same as for the left panel and shows which electrons (i.e. propagating along which open magnetic field 
lines) contribute to different parts of the spectrum. 
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Figure 2: Photon maps and light curves obtained for MSP with P = 5ms and B pc = 3.5 x 10 8 G. Top panel 
shows pulsar with inclination of 50°, while bottom panel shows the case for 70°. In the photon maps a 
photon density is colour coded in arbitrary units. On the x-axis a pulsar rotation phase <j> is presented while 
the y-axis shows a viewing angle £ in degrees. Pulse profiles for observers with viewing angles £ = 50° and 
70° are presented. 
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